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Fxeculivc  Sum  maty' 

The  work  i>erforimd  in  the  reporting  period  has  been  focused  on  continuation  of  Task 
l.l.  and  initial  work  on  Task  2.2  and  Ta.sk  2.4,  as  described  in  the  project  proposal.  The 
activities  related  to  Task  l.l  are  construction  of  a  representative  volume  element  (RVE) 
containing  disordered  fiber  distributions,  simulating  manufacturing  effects,  and  its 
computational  micromechanics  failure  analysis,  Ta.sk  2.2  is  concerned  with  RVE  level 
modeling  by  peridynamics  and  Task  2.4  addresses  fatigue  damage  with  peridynamics. 
Initial  results  for  all  tasks  are  reported  and  the  ongoing  work  is  outlined. 
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Task  1,1  Micro- level  crack  Initiation 


Background  and  motivation 

Manufacturing  of  polymer  matrix  composites  (PMCs)  results  in  disordered  fiber 
distributions  and  defects  such  as  voids  in  the  matrix.  We  focus  in  this  task  on  the  local 
stress  states  in  the  matrix  that  become  triaxial  under  any  remotely  applied  stresses 
resulting  from  the  service  environment  in  which  a  given  composite  structure  operates. 
The  local  stress  states  are  responsible  for  the  precursor  mechanisms  that  initiate  cracks. 
This  task  studies  the  point-failure  processes  that  become  critical  under  favorable  energy 
conditions.  The  two  basic  energy-driven  processes  are  cavitation  and  shear-band 
formation  in  polymers.  The  former  is  governed  by  the  dilatation  energy  density,  while  the 
latter  requires  energy  of  distortion  for  its  initiation. 

The  work  conducted  in  the  reporting  period  has  focused  on  cavitation  resulting  in  regions 
of  fiber  clusters,  where  favorable  conditions  for  cavitation  are  likely.  The  remote  loading 
considered  in  this  work  is  tension  normal  to  fibers  (so-called  transverse  tension). 

Approach  and  Results 

The  approach  taken  here  is  to  consider  the  type  of  fiber  clu.sters  that  arc  likely  to  result 
from  the  manufacturing  process.  In  most  cases,  the  resin  llovv  displaces  fibers  that  are 
initially  relatively  ordered  and  are  in  contact  with  one  another.  The  displacement  caused 
by  the  resin  How  induced  forces  will  be  the  outward  radial  movement  and  the  in-plane 
rotation  of  fibers.  Both  these  displacement  variables  will  be  random,  as  the 
manufacturing  process  cannot  be  controlled  in  a  deterministic  way.  In  the  simulation,  a 
pre-spec i lied  range  in  both  variables  is  implemented  by  using  uniform  random 
distributions. 

Fig.  1  (a)  below  shows  an  initial  (dry  )  cross-section  of  a  circular  region  of  a  fiber  bundle, 
and  Fig.  I  (b)  displays  the  disordered  fiber  distribution  after  the  random  displacements 
have  been  enforced.  As  can  be  seen,  the  post- manufacturing  fiber  distribution  has 
clustered  and  diffused  sub-regions.  As  reported  in  the  P*  Quarterly  Progress  Report,  the 
constraint  caused  by  closeness  of  the  fibers  results  in  favorable  conditions  for  cavitation. 
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(a) 


(b) 


Fig.  I .  (a)  Initial  region  of  a  dry  fiber  bundle  where  all  fibers  are  in  contact,  and  (b) 
Simulated  region  of  disordered  fiber  distribution  using  random  outward  radial  movement 
and  in-plane  rotation  of  fibers. 

Several  realizations,  such  as  that  shown  in  Fig.  I  (b),  were  generated.  The  number  of 
fibers  used  in  each  case  was  91.  Each  RVE  realization  was  analyzed  with  a  finite  element 
model,  where  the  local  stress  states  in  the  matrix  were  calculated  and  were  utilized  in  the 
failure  analysis. 

The  failure  analysis  consisted  of  determining  which  of  the  two  mechanisms  (cavitation 
and  shear-banding)  are  likely  to  initiate  first.  Figure  2  shows  the  dilatational  and 
distortional  energy  density  variations  with  the  remote  tensile  strain.  The  zero  strain 
corresponds  to  thermal  cooldown. 
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Figure  2.  The  variations  with  applied  strain  of  the  dilatational  energy  density  (lower 
curve)  and  of  the  distortional  energy  density  (upper  curve).  The  criticality  is  reached  at 
approximately  0.7%  strain  when  the  dilatational  energy  density  reaches  the 
experimentally  determined  critical  value  (0.2  MPa). 
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To  validate  whether  the  critical  dilatatlonal  energy  density  represents  hydrostatic  tension, 
the  three  principal  stresses  are  calculated  and  their  pair-wise  ratios  are  plotted  in  Fig.  3. 
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Figure  3.  The  pair-wise  ratios  of  principal  stresses  and  the  ratio  of  the  mean  to  max 
principal  stress  corresponding  to  the  results  in  Fig.  2  are  plotted.  At  0,7%  strain  the 
triaxial  stress  state  is  not  near  enough  to  a  hydrostatic  state  to  favor  cavitation. 


The  results  showm  in  Fig.2  and  Fig.  3  are  for  one  realization  of  the  RVE  simulating  a 
manufacturing  induced  disorder  in  fiber  distribution.  This  realization  illustrates 
unfavorable  conditions  for  cavitation,  A  total  of  six  realizations  were  performed.  Two  of 
the  six  did  fall  in  a  favorable  region  of  cavitation. 


The  effect  of  cooldown  temperature  from  processing  to  room  temperature  was  also 
studied.  The  results  above  are  for  82*^0  cooldown  temperature,  which  is  typical  for  glass- 
epoxy  composites.  However,  this  temperature  can  be  varied.  Today,  resins  are  available 
to  do  room  temperature  curing  of  thermosetting  resins.  Therefore,  the  cooldown 
temperature  was  varied  to  see  its  effect  on  initiation  of  cavitation  (and  consequently, 
fiber-matrix  debonding).  Following  summarizes  the  computed  results. 

When  the  cooldown  temperature  was  reduced  to  75*^C,  the  favorable  conditions  for 
cavitation  were  found  to  occur  earlier,  at  0.6%  strain.  Further  reducing  the  cooldown 
temperature,  however,  did  not  favor  cavitation.  On  the  other  hand,  increasing  the 
cooldown  temperature  required  increasing  the  mechanical  load  for  cavitation  to  initiate. 
These  results  have  important  implications  on  the  effect  of  manufacturing  on  failure 
initiation  under  transverse  loading. 
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Ongoing  work 

The  effect  of  manufacturing  defects  on  crack  initiation  will  continue  a  little  into  the  third 
quarter.  Then,  the  focus  will  be  on  transverse  crack  formation  and  its  multiplication  in  the 
presence  of  manufacturing  defects. 

Task  2.2  RYE  level  modeling  with  peridynamics,  and  Task  2.4:  Fatigue  damage 
modeling  with  peridynamics. 

Motivation 

Pores,  manufacturing  defects,  holes  and  round  notches  are  locations  where  fatigue  cracks 
may  initiate  from  or  arrest,  depending  on  the  loading  conditions.  Investigating  the 
evolution  of  fatigue  cracks  after  sinking  into  a  hole,  for  example,  is  critical  for  predicting 
damage  evolution  in  composite  materials,  in  which  pores,  inhomogeneities  and 
manufacturing  defects  play  critical  roles. 

Approach  and  Results 

We  have  introduced  a  fatigue  crack  growth  model  based  on  peridynamics  and  validated 
our  implementation  against  experimental  results  in  terms  of  the  crack  path  shape,  We 
perform  convergence  studies  in  terms  of  the  nonlocal  region  size  for  a  modified  compact 
tension  test,  which  leads  to  curved  fatigue  crack  paths,  The  experimental  crack  paths  arc 
strongly  inlluenced  by  the  location  of  a  hole  in  the  specimen,  and  the  peridynamic  results 
capture  this  sensitivity  well  (see  figure  4).  Fatigue  lifetimes  obtained  by  the  peridynamic 
model  are  in  good  agreement  with  experiments  for  different  crack  growth  rates  in 
different  cycle  ranges.  Different  from  methods  based  on  classical  continuum  mechanics, 
the  peridynamic  fatigue  crack  model  does  not  require  additional  criteria  to  guide  crack  or 
damage  growth.  In  particular,  we  use  the  peridynamic  fatigue  model,  w'ithout  any 
modillcations,  to  simulate  fatigue  crack  grow'th  in  a  two-phase  composite  in  which 
several  crack  initiation  points  exist  and  where  fatigue  crack  paths  interact  in  complex 
ways.  We  have  published  these  results  in  reference  1 1  ]. 
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Fig.  4.  Crack  paths  for  the  different  modified  CT  specimens  obtained  by  peridynamic 
model.  The  difference  between  the  samples  is  the  location  of  the  added  hole  (see  upper 
figure).  Results  published  in  [  1  ]. 

To  further  test  the  model  for  fatigue  crack  growth  after  a  crack  sinks  into  a  hole,  we 
proposed  improvements  that  add  the  fatigue  limit  not  only  to  the  initiation  phase,  but  also 
to  the  propagation  phase.  We  demonstrate  that  the  model  simulates  all  three  phases  of 
fatigue  failure  (initiation,  propagation,  and  final  failure)  with  an  example  in  which  a 
fatigue  crack  sinks  into  a  cutout  and  re-initiates  from  a  different  location  along  the  cutout, 
grows,  and  lead  to  final  failure  of  the  structure  (see  Figure  5).  The  fatigue  crack  path 
from  the  improved  model  agrees  with  an  analysis  based  on  strain  concentrations. 
Convergence  studies  show  that  the  peridynamic  results  are  correct  once  the  nonlocal  size 
is  smaller  compared  with  the  size  of  relevant  geometrical  features.  These  results  will 
appear  in  publication  [2]. 
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Fig.  S.  Crack  nucleation  after  the  original  crack  sinks  into  the  added  hole.  Damage  index 
map  showing  the  transition  between  fatigue  crack  growth  and  quasi-static  fracture/final 
sample  failure.  Results  to  appear  in  [2]. 


Applied  to  a  fiber-reinforced  polymer  composite  material,  the  fatigue  crack  model 
produces  results  consistent  with  those  seen  in  experiments  (see  Figs.  6  and  7).  In  Fig.  6, 
we  show  a  sample  RVE  in  which  the  red  colored  disks  are  the  fiber  cross-sections  and  the 
blue  is  the  polymer  matrix  material. 


Fig.  6.  Cross-section  of  an  RVE  in  an  FRC.  Fatigue  loadings  are  applied  in  the  vertical 
direction. 
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The  results  in  Fig.  7  show  that  fatigue  cracks  can  initiate  at  multiple  sites.  The  local 
arrangement  of  fibers  leads  to  stress  concentrations  that  produce  matrix  micro-cracks  that 
join  together  in  generating  matrix  macrocracks.  Our  peridynamic  model  captures  this 
very  well.  It  is  interesting  to  also  observe  that,  in  certain  cases,  cracks  grow  around  the 
fibers,  indicative  of  fiber-matrix  debonding  failure.  In  the  simulations  shown  in  Fig.  7, 
the  interface  is  slightly  “stronger'  than  the  matrix,  and  because  of  that,  debonding  is  less 
dominant  here.  Our  further  studies  will  focus  on  the  interface  fatigue  properties. 


Fig.  7.  Damage  maps  at  different  fatigue  cycles  showing  fatigue  crack  nucleation  at 
multiple  sites,  growth  and  coalescence  of  cracks.  Fiber-matrix  debonding  is  also 
observed. 
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